Ancilla-less selective and efficient quantum process tomography 
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Several methods, known as Quantum Process Tomography, are available to characterize the evo- 
lution of quantum systems, a task of crucial importance. However, their complexity dramatically 
increases with the size of the system. Here we present the theory describing a new type of method 
for quantum process tomography. We describe an algorithm that can be used to selectively estimate 
any parameter characterizing a quantum process. Unlike any of its predecessors this new quantum 
tomographer combines two main virtues: it requires investing a number of physical resources scal- 
ing polynomially with the number of qubits and at the same time it does not require any ancillary 
resources. We present the results of the first photonic implementation of this quantum device, char- 
acterizing quantum processes affecting two qubits encoded in heralded single photons. Even for this 
small system our method displays clear advantages over the other existing ones. 
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Introduction. Quantum process tomography (QPT) 
[MTU] is a task of fundamental and practical importance. 
In fact, it can be used not only to study the evolution 
of generic quantum systems but also provides the nec- 
essary information to characterize the most important 
noise sources affecting a quantum information proces- 
sor. This task is essential to design good error correction 
strategies which must be used in order to prevent the de- 
coherence induced by the coupling with the environment. 

The full characterization of a quantum process (i.e., 
full QPT) is always a hard task that requires invest- 
ing resources that scale exponentially with the number 
of qubits (such number is hereafter denoted as n, being 
D = 2 n the dimension of the corresponding set of states). 
As the complete characterization of a quantum process 
requires 0(D 4 ) real numbers (see below), it is necessary 
to have methods that enable a partial characterization 
of quantum processes in an efficient manner (i.e., invest- 
ing resources scaling polynomially with the number of 
qubits of the system). Ideally, two requirements must be 
imposed on a good method for QPT: First, the method 
should enable us to select only a few parameters of the 
process to be efficiently determined. Also, the method 
should avoid the use of expensive resources such as clean 
qubits (i.e., qubits that are required to be perfectly iso- 
lated from the environment even though they interact 
with the system) . Existing methods for QPT do not sat- 
isfy these criteria: They are either inefficient or they use 
expensive ancillary qubits. In this report we present the 
first method satisfying the two requirements. In fact, our 
method can be used to perform partial QPT without us- 
ing any ancillary qubits while efficiently estimating any 
parameter of the \ matrix of a quantum process (with 
accuracy independent on the number of qubits of the sys- 
tem). 

In essence the method maps any relevant parameter of 
a quantum process onto the average transition probabil- 



ities between a special set of quantum states through a 
quantum process. Such states can be efficiently prepared, 
sampled and detected. Using a heralded single photon 
source we fully implement this method to perform QPT 
on any process jointly affecting the polarization and mo- 
mentum qubits of single photons. 

Ancilla-less quantum process tomography. An un- 
known quantum process affecting a physical system com- 
posed of n qubits can be represented as a linear map 
(a quantum channel) taking initial states po into final 
states p = A(po)- Any such map can be parametrized 
in terms of a x _m atrix as follows: Given a basis of op- 
erators E a (with a — 0, ...,D 2 — 1) any channel can be 
written as A(p) = ^ b XabE a pEl The D 2 x D 2 matrix 
X parametrizes the map which. It is hermitian and trace 
preserving if and only if \ = X f and Y, a b XabE ] h E a = /. 
Positivity of the x _m atrix is equivalent to the complete 
positivity (CP) of the map A. Simple counting arguments 
show that the number of real parameters defining such 
map is 0(D 4z ). This makes the task of full QPT hope- 
less for large systems (unattainable even for systems of 
a few tens of qubits). Moreover, until recently, methods 
available for QPT [1 were such that the evaluation of 
a sub-set of elements Xab also required resources scaling 
exponentially with n. In fact, the standard approach to 
QPT consists of estimating the transition probabilities 
Pkk' = Tr(pk'A(pk)) for a complete set of initial and fi- 
nal states pk and pw (with /c, h' = 0, ...,D 2 — 1). After 
this exponentially large number of experiments the Xab 
elements are determined by inverting a set of linear equa- 
tions which is also exponentially large. Improvements 
in this inefficient method were presented recently [3H5]. 
Such methods enable the efficient estimation of any di- 
agonal element Xaa- 

However, to estimate off-diagonal 
elements Xab such methods require the use of extra an- 
cillary clean qubits. We now present a general method 
that overcomes this problem. It is efficient and it does 
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S = {\<f>j),j = 1, is a 2-design 



FIG. 1: Experimental setup. Single 810nm heralded pho- 
tons are generated at a BBO cristal by Type II parametric 
down conversion. Arbitrary states are prepared and measured 
with a combination of single qubit unitary gates and a con- 
trolled operation as exemplified in the inset, colours show the 
correspondence between circuit and physical implementation. 
Path-qubit unitary gates are controlled with phase plates and 
Sagnac interferometers. Polarization-qubit unitary gates are 
achieved with a combination of wave plates affecting both 
paths. A controlled operation with control on path and target 
in polarization is implemented with wave plates at different 
angles in each path. The different processes studied were set 
up in the zone marked as S. Single-Mode fibers clean the 
photon's spatial mode ensuring good interferometer visibility 
(> 95% in the Sagnacs and w 90% in the Mach-Zehnder). 
Not shown in the figure is the Mach-Zehnder 's active phase- 
sensitive stabilization mechanism (see supplementary mate- 
rial). 



not require the use of extra clean qubits. The method is 
based on an important fact: any element of the ^-matrix 
can be interpreted as the average survival probability (i.e. 
fidelity) of a certain quantum map [4]: 
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This is the average over the entire Hilbert space of the 
survival probability of a map K a b defined as K a b{p) = 
k(E\pEb) (i.e., it is obtained by first transforming p into 
E\pEb and then applying the channel A). The efficient 
estimation of F a ij is equivalent to that of Xab- How- 
ever, two main obstacles are apparent impediments for 
the efficient estimation of F a b- The first is that aver- 
aging over the entire Hilbert space apparently requires 
preparing and measuring a infinite number of quantum 
states. The second obstacle is that the effective channel 
A a 5 is not physical (it is generally not a CP map un- 
less E a = Eb). The first obstacle can be surmounted by 
using the tools presented in [3j [4] . Thus, we can trans- 
form the integral over the entire Hilbert space into a sum 
over a finite set of states that form a so-called 2-design, 
which exist for any dimension [TTHT6] . In fact, if the set 
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The exact computation of Xab involves finite but expo- 
nentially large resources since K = 0(D 2 ). However, by 
randomly sampling over a subset of the 2-design, after M 
experiments one estimates F a 5 with an error that scales 
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M \ 



1 



M-l 
K-l 



. The error scales roughly as 



1/VM for M < K and vanishes for M = K [3H5]. Thus, 
the precision fixes the required number of experiments, 
not the size of the Hilbert space. 

The way to surpass the second obstacle, i.e. the fact 
that A a 5 is not a physical map, is to notice that it can 
be obtained as the difference between two CP maps. To 
be precise, let us describe how to achieve the efficient 
estimation of the real part of Xab- We exploit the con- 
nection with the real part of F a b shown in (2) and define 
the fidelities F^ b of two efficiently obtainable CP maps 
(see below) as: 

Fab = £<^|A((£a ± E^\<P 3 ){^\(E a ± E b ))\<p 3 ). (3) 



The desired fidelity is obtained by measuring F^ b and us- 
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F„ b . Therefore the estimation 



ing that 2Re(F ab ) 

of F a 5 (and with it, the estimation of Xab) 1S summa- 
rized as follows: 1) Randomly choose an element of the 
2-design \4>j)] 2) Efficiently prepare the state obtained 
by acting with (E a ±Eb)^ on the state chosen in the first 
step (see below); 3) Apply the channel A to the resulting 
state; 4) Estimate the probability to detect \<pj) as the 
output state. By repeating this process M times we esti- 
mate F^ h with an accuracy scaling as 1/ y/M. This is the 
core of the method that is the first one enabling the effi- 
cient estimation of any x~ m atrix element without using 
extra ancillary resources. The above steps are applica- 
ble because: i) the states of the 2-design, \4>j), can be 
efficiently prepared and detected; and ii) the states ob- 
tained by acting with (E a ±Ei))^ on \<pj) can be prepared. 
The details of this preparation process and other techni- 
cal remarks including determination of the imaginary x 
elements can be found in the supplementary material. 

We implemented this method in an experiment to fully 
characterize several quantum channels affecting n = 2 
qubits. As a 2-design we used the D{D + 1) = 20 eigen- 
states of D + 1 = 5 mutually unbiased bases (MUBs). 
In particular we chose the three separable bases whose 
generators are X, Y and Z for each qubit and two en- 
tangled ones generated by the operators {X <g>Y,Y ® Z} 
and {Y ®X,Z®Y}. 

Photonic implementation. Several methods for com- 
plete or partial quantum process characterization have 
been demonstrated in different experimental setups [2, 
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FIG. 2: Full reconstruction of two channels. The x matrices for a) the identity process and b) a controlled U c process 
(see text). Shown are the real and imaginary parts of each one all in excellent agreement with the expected theoretical results. 
The fidelity between these reconstructed processes and the ones reconstructed by the standard inefficient Nielsen and Chuang 
method are 91% and 93% respectively. 



[171420] . Here we implement our QPT method, which 
can selectively determine any parameter characterizing a 
quantum process, using a heralded single photon source 
obtained by Type-II parametric down conversion at a 
BBO crystal with a 405nm diode laser producing twin 
photons at 810nra. Single photons encode two qubits: 
one associated to the polarization degree of freedom and 
the other to the path. As described in Figure [I] the exper- 
iment is divided in three stages: i) state preparation, ii) 
evolution with the quantum channel to be tomographed 
and iii) transition probability measurement. Polarization 
qubits are controlled with several wave plates while path 
qubits are controlled by three interferometers. Controlled 
operations are done with wave plates at different angles 
in each path. By appropriately combining a self-stable 
Sagnac interferometer, phase plates and wave plates we 
can prepare any desired state and measure any of the 
states of the 2-design as required [21 . Figure [I] and its 
inset show the circuit equivalence of each part. State 
preparation and measurement are connected by an ac- 
tively stabilized Mach-Zehnder interferometer in which 
the process is embedded (see supplementary material). 

Various processes were studied: the identity; a unitary 
on the polarization qubit (a wave plate on both paths); 
an operation U c in which a different unitary is applied 
to the polarization qubit depending on the path qubit 
(U c = (I - Z) <8> Z/2 + (I + Z) <8> X/2 is implemented 
with wave plates at different angles in each path) and 
a noisy version of U c (noise in the path qubit is added 
by sweeping the Mach-Zehnder' s phase). Figure [2] shows 
the full reconstruction of the x matrix for the identity 
and U c . In all measured processes we obtained excellent 
agreement with x~matrix of the ideal process and with 
the one measured by the standard method pQ. We com- 
puted the fidelity between the x matrix obtained with 
our method and the one measured using Nielsen and 



Chuang's method. We obtained that for all the imple- 
mented channels such fidelity is above 90%. For full 
QPT we independently measured the values of all Xab 
elements. To do so, for each element we prepared the 
20 states in the 2-design and measured their survival 
(and non-survival) probabilities. A simple algorithm to 
prepare states (E a ± Eb)\(/)j) was developed. Full char- 
acterization of the 256 elements of the x matrix involves 
256 x 40 = 10240 transition probabilities. Fortunately 
many of them coincide and the number of different tran- 
sition probabilities is much lower. In our case, a full char- 
acterization of a two qubit process requires only 140 dif- 
ferent experimental settings, each giving 4 probabilities, 
resulting in 560 transition probabilities to be measured 
(see supplementary material). 

However for the above full QPT, we do not take ad- 
vantage of the most powerful aspect of this method: effi- 
ciency and selectivity. We also performed efficient partial 
quantum process tomography measuring useful proper- 
ties of the channel without fully determining the x ma- 
trix. In this way we clearly show the advantage of our 
method over previous ones. Suppose we are interested 
in determining how close a given process is to a tar- 
get process. A good measure of such distance is pro- 
vided by the average fidelity of the channel obtained by 
composing the inverse target operation and the chan- 
nel A. It is simple to show that this average fidelity is 
F = (DTr(xx) + ])/iP + 1)? where x is the x _ matrix 
of the inverse target channel. Such matrix can be simply 
obtained analytically and is typically small: For example, 
for the identity channel the only non- vanishing element 
is Xoo = 1- Therefore to estimate the fidelity of the iden- 
tity we only need to estimate xoo- On the other hand, 
for the above mentioned controlled operation U c , the x 
matrix has 16 non-vanishing elements (4 diagonal and 
12 off-diagonal ones). Therefore, to measure the fidelity 
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I | Identity Q Controlled U c gate Q Controlled U c gate with dephasing on one qubit 



FIG. 3: Efficient parameter estimation. By determining 
only some parameters characterizing a quantum channel dif- 
ferent relevant questions can be answered. These parameters 
can be determined exactly or approximated by sampling over 
only some of the possible measurements needed to fully de- 
termine them. Depending on the complexity of the question 
more or less measurements are needed. In any case these val- 
ues can be determined statistically with an error scaling with 
the amount of measurements better than 1/y/M eliminating 
the need for exponential amount of measurements. Shown are 
the fidelities of the measured processes to two different target 
processes: the identity (left) and a controlled gate U c (right) 
as a function of the amount of measurements done and for 
ten different choice of sampling order. Results for different 
channels show how quickly the estimation can differentiate 
between each one and converge to a value close to the exact 
one. Also shown in black and coloured dashed curves is the 
statistical maximum deviation expected for each result. 

of such channel we only need to estimate 16 elements 
of the x _ma trix. Moreover, such coefficients can be es- 
timated with increasing precision by increasing the size 
of the sample. In Figure [3] we show how these fidelities 
converge when the sample size is increased (curves corre- 
spond to different random choices for the order in which 
we sample over the 2-design). In all cases we see that 
it is possible to decide if the channel A is close enough 
to the target channel by making a number of measure- 
ments that is much smaller than the one required for full 
QPT. In such figure we also show that the same method 
reveals the presence of noise in the controlled operation. 
Not only an exact answer to such questions requires fewer 
resources than in previous methods [TJ [6] but a good es- 
timate can be obtained with a number of measurements 
that does not scale exponentially. 

Conclusion. The method presented accomplishes one 
of the fundamental tasks needed to build a general pur- 
pose quantum computer: it can efficiently and selectively 
characterize any quantum process. The power of this 
quantum device arises from the ability to prepare, sam- 
ple and measure a special set of states. We showed that 
even for a two qubit system the strategy we presented 
displays clear advantages over preexisting methods. Our 
photonic experiment confirms this by identifying crucial 
information determining the nature of different quantum 
processes in a selective manner. The crucial task of esti- 
mating how close a certain evolution resembles a target 



process can be performed with our technique in a selec- 
tive and efficient way. This, together with the estimation 
of the most important sources of errors affecting a quan- 
tum process are examples of the type of usage that this 
type of quantum tomographer will have. 
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Suplementary Material 

Calculation of probabilities. For each configuration, 
clicks in each output were recorded during 10 seconds. 
Probabilities were calculated by normalizing this amount 
of clicks to a reference measured with the identity process 
and by preparing and measuring in the computational 
(ZZ) base. This reference was updated in periods of less 
than 25 minutes to account for slow fluctuations in the 
laser power. All optical fibers on the signal side where 
single mode fibers. The average amount of maximum 
counts per second were: 1800000 on the herald, 2600 on 
the signal and 600 coincidence counts. 

Interferometer Stabilization. The Mach-Zehnder in- 
terferometer was actively stabilized by a phase-sensitive 
closed-loop system. An auxiliary 405nm laser runs par- 
allel to the 810nm single photon beam path and its in- 
terference signal is m easured with a photo-diode. Using 
this signal the position of one of the mirrors in the inter- 
ferometer is actuated with a piezoelectric disc. The ro- 
tation of the 810nm wave plates only affects the 405nm 
reference signal in that its intensity varies considerably, 
but not the position of the maxima and minima. By 
making the piezoelectric disc oscillate with a small am- 
plitude (~ lOnra) at a frequency of ~ 3kHz and com- 
paring the relative phase of the measured oscillation with 
the forcing one, an analog circuit determines the neces- 
sary correction to the offset of the oscillation so as to 
maintain the interferometer always at a minimum of the 
interference fringes of the 405nm reference beam. With 
this simple phase sensitive technique we get stability of 
approximately A/20 for the single photons. 

Imaginary parts of Xab • In order to obtain the imagi- 
nary part of an off-diagonal element, we should consider a 
slight modification of the real part scheme shown above. 
First, we must consider the following two CP channels: 

3 

Then it is strightforward to obtain the imaginary part 
of one of the fidelities from equation 1 by measuring F^ h 
and considering 2Im(F a b) = F^ b — F~ b . 

Efficient state preparation. It is important to use a 
special 2-design adapted to the basis of operators E a . 
We choose E a as generalized Pauli operators built as n- 
fold tensor product of the identity I or one of the three 
Pauli operators (X, Y or Z) on each qubit. These oper- 
ators form the Pauli group that has D 2 elements (up to 



phases) . This group can be partitioned into D + 1 com- 
muting subgroups each of which contains D operators 
(including the identity) which are obtained as all possi- 
ble products between n independent generators. Each 
commuting subgroup defines an orthonormal basis of the 
Hilbert space, formed by the eigenstates of the opera- 
tors in the set. These (D + 1) bases are mutually unbi- 
ased (MUBs) [22} [23] : any state of one basis is an equally 
weighted superposition of all the D states of any other 
basis. The set of all D(D + 1) states belonging to the 
(D + 1) bases form a 2-design with special properties. 
This is the 2-design we will use in our method, and we 
will denote it as S = {|0^),a = 0,...,D,i = 1, ...,£>}. 
The index a labels the different MUBs and the index 
i labels each state in each basis. For our method we 
will use the following properties of this 2-design: a) Any 
state \4 a) ) 

can be generated from any computational 
state (i.e. a joint eigenstate of all Z{ operators) by an 
efficient quantum circuit [4]. b) Any Pauli E a is such 
that E a \(j)f^) = |<A^), that is to say, Paulis are trans- 
lations within each basis. The transition rule, i.e., the 
expression that determines i' as a function of (z,a,a) 
can also be efficiently obtained. In fact, such expression 
depends only on the commutation (or anti-commutation) 
relations between E a and the n operators that are chosen 
as generators of the basis a. c) Also the normalized state 

I*±?o,m) = K (E a ±E h )\(j) { f ) ) can be generated efficiently 
from any computational state. The simplest way to do 
that is to prepare first a superposition of appropriately 
chosen computational states and later apply a change of 
basis. These tasks can be efficiently performed. The nor- 
malization constant K is also efficiently computable. 
We need to prepare states of the form (E a + e % PEb)\(jff^), 
where f3 is a multiple of tt/2 (odd multiples are required 
for the measurement of the imaginary part of \ab) an d 
|</>^) is one of the states from the 2-design. To do this, 
we first fix an ordering of the states within each basis. 
On the computational basis (a = 0), we choose the lex- 
icographic ordering. For any other basis we will use the 
convention \<pf) = V$\(j)f), where Vq is the correspond- 
ing change of basis operator from |4 . The states we 
prepare are then of the form 

(Ea + e^E^VfX^l^) (5) 

where |0o°^) is the vector of the computational basis that 
has all zeros and is an operator that has X on each 
qubit where the binary decomposition of i has a one. 

Since V$ is a Clifford group operator built with 0(n 2 ) 
Hadamard, CNOT and phase gates, it is efficient to com- 
pute how E a and Eb transform into E a and Eb under 
conjugation via ^[3 . This yields: 

^(ijW+eW))^) (6) 
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And since the application of a Pauli operator on a com- 
putational basis state yields another state from that ba- 
sis, the required state can be restated as 

K?(lO+e* 7 |</>°» (7) 

which is the change of basis circuit acting on a state that 
is efficiently prepared via a Hadamard gate, 0{n) CNOT 
gates and at most three phase gates. The normalization 
constant is readily obtained from Q as the norm of the 
state prior the application of the change of basis. 

Only 560 probabilities. It is no surprise that full pro- 
cess tomography requires an exponential number of prob- 
ability measurements. Just for diagonal tomography, 
each single coefficient is an average of 20 probabilities. 
Since there are 16 such coefficients this would require 



on the order of 320 probability measurements for full di- 
agonal tomography. However, many of these probabili- 
ties are repeated. For instance, if we were to measure 
the diagonal x coefficients corresponding to the opera- 
tors X (g) X and Y (g) Y, it is straightforward to see that 
both operators acting on the state |0q ^ ) yield the same 
state, up to a phase. When it comes to off-diagonal to- 
mography, many more of those probabilities are repeated. 
Since the \ matrix is defined by 0(D 4 ) real numbers and 
each each of those requires 0(D 2 ) probabilities to be ob- 
tained, it can be seen that full process tomography will 
require between 0(D 4 ) and 0(D 6 ) probability measure- 
ments, both exponential on the number of qubits. In the 
2 qubits case it was found out that measuring only 560 
probabilities was enough for full process tomography. 



